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One of the more intriguing hallmarks of many-body 
systems is that at zero temperature quantum fluctua- 
tions can drive the system to a drastic change of state, 
commonly known as a quantum phase transition (QPT). 
A paradigmatic model for QPTs is the one-dimensional 
Ising model Recently, experimental realizations of 
one-dimensional spin chains have been suggested, where 
a quantum simulation of the system close to the phase 
transition is possible, and a wide freedom on the con- 
trol of the parameters is achieved 043- The quantum 
control of many-body systems by a driving field has at- 
tracted considerable interest both theoretical and experi- 
mental with workers from very different communities be- 
ginning to look at driven models 0-[Il- The possibility 
to manipulate the quantum state of a system by means 
of a classical external control allows one to explore novel 
states of matter and effective interactions which are ab- 
sent in equilibrium [l5l - [l8} . In the presence of an external 
control, quantum resonances and symmetries play an im- 
portant role (l9l - l22| . In particular, as a consequence of 
a generalized parity in the extended Hilbert space fl9| . 
under the effect of periodic driving the tunneling can be 
slow down or totally supressed in a perfect coherent way, 
a phenomenon commonly referred to as coherent destruc- 
tion of tunneling (CDT) [23|,|24(. Rather recently, the ex- 
tension of this concept to many-body systems has been 
addressed in the context of the Mott-superfluid transi- 
tion in ultracold s yste ms either theoretically [l2| as well 
as experimentally [lj] , and in a two- mode Bose-Hubbard 
model with time-dependent self- interaction strength [Io| . 

The dynamics of one-dimensional spin chains has been 
addressed extensively when the system is driven slowly 
through the critical point [25l - |27| , where there is a diverg- 
ing relaxation time and correlation length, and the dy- 
namics cannot be adiabatic in the thermodynamic limit. 
As a consequence of this, the final state of the system con- 
sists of ordered domains whose finite size depend upon 
the velocity of the parameter variation [28j . A nontrivial 
oscillation of the magnetizatio n | 29| and the connection 
between symmetry and CDT [30| has been investigated 
in a finite size periodically-driven Ising model. Further- 
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more, under the effect of a nonadiabatic external control 
of the transverse field, the Isin g ch ain exhibits dynami- 
cal freezing of the response (3fl. |32||. and synchronization 
with the external driving in the asymptotic dynamics as 
a consequence of destructive interference in time (33|. 

Our aim in this paper is to describe the nonequilib- 
rium behaviour of a one-dimensional Ising model under 
the effect of a nonadiabatic monochromatic transverse 
field from the perspective of quantum criticality. In par- 
ticular, we describe the dynamics by means of an effec- 
tive Hamiltonian which simulates the dynamics of an un- 
driven system. We show that in the asymptotic dynamics 
the nonequilibrium quantum phases correspond to states 
of the system which are synchronized with the driving. In 
contrast to previous works j3ll - l33| . however, we describe 
the role of many-body CDT in the critical behaviour by 
investigating signatures of criticality both in the labora- 
tory frame as well as in the rotating frame. 

The paper is organized as follows: In sec. fl] we dis- 
cuss the equations of motion and describe quantum res- 
onances by considering the excitation spectrum of the 
undriven system. In sec. |IT] we describe the physics of 
the system by means of the rotating wave approxima- 
tion and discuss signatures of criticality based on the 
description of the quasienergy dispersion and CDT. In 
sect. IIIII we describe the quantum dynamics in the lab- 
oratory frame by considering the quantum evolution of 
the system when is initialized in a paramagnetic ground 
state. Furthermore, we study signatures of criticality by 
considering cycle-averaged expectation values of physical 
observables. Finally, discussion of the results is presented 
in sec. IIVI 

I. QUANTUM RESONANCES IN THE ISING 
MODEL. 

The periodically-driven transverse Ising model de- 
scribes the dynamics of N interacting two-level systems 
in a time-dependent transverse local field [3ll - l33j 

N N 
i=l i=l 

where af are Pauli matrices at the ith site and we 
assume periodic boundary conditions af = cr'^ +1 for 



2 



a E {x,y,z}. In the following we shall consider J > 0, 
and a monochromatic modulation of the transverse field 
with a static contribution g(t) = go + gi cos fit. Via 
a Kramers- Wannier self-duality transformation [HI, HH , 
and after a spin rotation, it is possible to map Hamilto- 
nian Eq. ([TJ into a dual periodically-driven Ising model 

N N 

H D (t) = -^)^rfrf +1 -J^rf, (2) 

i=l i=l 

where T 4 z rf +1 = erf, and crf^af = rf. Therefore, the 
study of the nonequilibrium QPT in the Hamiltonian 
Eq. (TTJ), allows one to get a physical picture of the cor- 
responding QPT in the dual model where the spin-spin 
interaction is time-dependent. 

Associated with the Hamiltonian Eq. (JlJ is a conserved 
parity ft, such that [H(t), ft] = 0, which is given by 

N 

ft = (g)af. (3) 

i=i 

In the thermodynamic limit N — > oo, the undriven 
Ising model (gi = 0) exhibits a second-order QPT at 
<7o = J from a symmetric paramagnetic phase (go > J) 
to a symmetry-broken ferromagnetic phase (go < J) [J- 
Our aim in this paper is to study the new aspects of 
criticality under the effect of driving. 

In this section we provide the basics on the formal- 
ism used to describe the time-dependent Ising model. In 
particular, we find a resonance condition related to m- 
photon processes under the effect of driving. 

A. The dynamic Bogoliubov de Gennes equations 

In this paper we consider the restriction of the Hamil- 
tonian Eq. ([I} to the subspace with even (+) num- 
ber of fermionic quasiparticles (see appendix A). After 
a Jordan- Wigner transformation and a discrete Fourier 
transform of Hamiltonian Eq. ([!} we obtain 

H(t) = [2 (g(t) - Jcosk) (c\c k + cl k c- k ) - 2g(t) 

k>0 

+ 5^2JsmA(4c_ fc + c_ fc Cfc) = 5^#fc(t) I (4) 

fe>0 k>0 

where c£ and c k are fermionic operators [l[ . For finite size 
N of the spin chain, the quasimomentum is restricted to 
k e {±%,±?fr, . . . ,± (JV ~ 1)7r }. In the following, we focus 
on the thermodynamic limit N — > 00, where we have 

— 7T < k < 7T. 

Even if we prepare the system initially in a ground 
state of the undriven model, under the effect of nonadi- 
abatic external driving, the system will experience tran- 
sitions to excited states. Rather recently, a formalism 
has been developed to deal with this kind of dynamical 



situations [l3|, 1251427} . The idea is based on the fact that 
Hamiltonian Eq. Q conserves momentum and parity. 
As a consequence, we can use the BCS ansatz for the 
evolution of the quantum state of the system 

|V I t) = 0(ufc(*)|l-fc,l fc ) + t;fe(t)|O_ fc ,O fc » J (5) 

fe>0 

which implies that for a given quasimomentum k, the 
quantum evolution is restricted to the Nambu sub- 
space {|l-fc, lfe), |0_k, Ok)}, consisting of doubly occupied 
l_fc, It,) and unoccupied 1 _ ^ , ^ ) states of ±k fermions 

niEim. 

The matrix representation of the operator Hk(t) in the 
Nambu subspace is given by 

TT (t\-( - 2oj k A fe \ 

Hk ®-{ A k -At))' (6) 

where H k (t) is the Bogoliubov de Gennes (BdG) Hamil- 
tonian, cjf. = 2Jcosfc, Aft = 2Jsinfc, and fx(t) = 2g(t). 
By defining the spinor ^l(f) = (u^(t), v^(t)), and con- 
sidering the BCS ansatz Eq ([5]), it is possible to show 
that the coefficients u k (t) and v k (t) should satisfy the 
differential equation 

i±* k {t)=H k Wk(t), (7) 

which constitutes the dynamical version of the Bogoli- 
ubov de Gennes equation (2f|. At this point we have 
translated the many-body problem into the solution of 
the time-dependent Schrodinger equation for an effective 
two-level system. 

Under periodic driving, the Floquet theorem states that 
the solution of Eq. ([7]) can be written as 

* fc (t) = A+e- ie i +h $ k + \t) +A_e- is i~ h $ k - ) (t), (8) 

where $1 (t) denote the Floquet modes corresponding 

to the quasienergies —0/2 < < Q/2. Furthermore, 
in the extended Hilbcrt space 1Z (E> T, where 1Z is the 
Hilbcrt space of square integrable functions and T the 
space of time-periodic functions, the Floquet states sat- 
isfy the eigenvalue problem 

% k ^\t)=e k X) ^\t), (9) 

where A e {+,-}, T~Lk = H k (t) - U k -§- t is the Floquet 

Bogoliubov de Gennes (FBdG) Hamiltonian, are 

the quasienergies, and the Floquet modes (t + T) = 

<5> k X) (t) have the same period T = 2tt/H, as the external 
driving [H-IH. 

B. Resonance conditions 

In the thermodynamic limit, the Ising model is char- 
acterized by an infinite set of collective excitations. Un- 
der the effect of an external driving, the possibility of 
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multiphotonic resonances arises 

Ql IH H2|. To study 
such quantum resonances, let us consider the system in 
the absence of driving g\ = 0. In this case, the Flo- 
quet modes and the quasienergies become the stationary 
states <jy^ and the excitation spectrum ejj^ = — uik ± £fc 
of the undriven system, respectively [UTS]. Therefore, 
the solution of Eq. ([7]) can be written in the form of Eq. 
(HI) as follows 



ip k (t) = a+e ie * 



The energy gap is given by AEk 
where 



a_e 



(10) 



- ,(+) 



c(-) - 



£fe = 2 a/ (go — Jcosk) +(Jsmky 



= 2e fc , 



(11) 



In the semiclassical theory of light-matter interaction, we 
can interpret a Floquet state as a light-matter quantum 
state containing a definite, though very large, number 
of photons [lj|. Multiple transitions between quantum 
states of the spin chain that are not directly coupled by 
the interaction can occur by means of intermediate states 
with a different number of photons present [l9l - [2l| . In 
particular, m-photon transitions occur when the condi- 
tion 



AE h = mil. 



(12) 



with integer m is satisfied. For a parametric oscillator 
with fundamental frequency Eq. (1121) is the usual 
resonance condition |3g]. In Floquet theory, Eq. ([12"]) 
implies the existence of a crossing between the single- 
particle energy levels e& and —ek when the energy spec- 
trum is folded into the Brillouin zone [19]. Such a cross- 
ing occurs at the wave vector 



± arccos 




(13) 



where the resonance condition is fulfilled, as is depicted 
in Fig. [T] (a). Fig. [T] (a) depicts the energy dispersion 
relation of the undriven system and the dashed lines in 
Fig. [U (b), the corresponding folding of the energy spec- 
trum into the first Brillouin zone —51/2 < < 51/2. 
In this paper we focus on the weak spin-spin coupling 
limit go,fl^>J. In this limit the multiphoton resonance 
condition reads 



So 



mil 



(14) 



Such resonance condition will be used in the next sec- 
tion to perform a description of the system based on an 
effective time- independent Hamiltonian which is valid for 
parameters close to a multiphotonic resonance. 



II. PHYSICS IN THE ROTATING FRAME 

The dynamical BdG equations allow us to investigate 
the full quantum evolution of the driven system. As we 
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FIG. 1: (Color online) (a) Typical gapped energy spectrum 
±efe of the undriven system corresponding to the paramag- 
netic phase go 3> J. In this case, the energy dispersion is 
slightly curved because the curvature is proportional to the 
spin-spin interaction strength J. (b) The dashed lines depict 
the spectrum ±tfc when it is folded into the first Brillouin 
zone, the crossing at k — ±fco is related to a 2-photon res- 
onance. The solid lines represent the quasienergy spectrum 
±£fe,m for a driving amplitude gi/Q = 1. The parameters are 
m = 2, J/Q = 0.01, and g /Q. = 0.505. 



are interested in the asymptotic quantum dynamics and 
the description of critical signatures, it is convenient to 
describe the dynamics of the system in a rotating frame. 
In the weak spin-spin interaction limit, it is possible to 
neglect the fast oscillations in the rotating frame, and 
a description of the system based on the description for 
time-independent systems is possible via an approximate 
effective Hamiltonian. 



The rotating wave approximation and the 
effective Hamiltonian 



Motivated by the m-photon resonance condition in the 
weak spin-spin coupling limit Eq. (|14[) , we perform a de- 
scription of the system based on the rotating wave ap- 
proximation (20| . Let us perform a unitary transforma- 
tion of Hamiltonian Eq. (jTJ) into a convenient rotating 
frame via the unitary operator 

U^\t) = exp (ia m (t) 5>? j = J] Ui m) (t) 

\ i=l / k>0 

= Y[ exp -2ia m {t)(c\c k + c'_ k C-k - 1) , (15) 
fe>o 
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where a m (t) = m(f2/4)t+|^ sin ilt. In the rotating frame 
the dynamics is governed by the Hamiltonian (t) — 
[UWQWHUWy), where H = = E fc >o#* = 



J2k>o(Hk{t) — Hk-§t) is the Floquet Hamiltonian. The 
explicit form of this operator is given by 



N j N T N 

£(m) (<) = _ 6 (m) J2af--{1 + cos[4a ro (t)]}£ of (7f +1 - - {1 - cos[4a m (t)]}£ 



! = 1 



J ^ J ^ 

+ - sin[4a m (i)] CT f CT f+i + 9 sin[4a m (i)] ^ af a.f +1 , 



(16) 



where the detuning = go — m(f2/4) describes how far 
the system is from resonance, and m is an integer that 
denotes the order of the multiphotonic resonance [Tg| . 
By using the identity 



exp(iz sin fit) 



E 



Ji{z) exp(iZQi), 



(17) 



where ^(z) is the Zth-order Bessel function [331, the 
Hamiltonian Eq. (|16p can be written in the form 



(18) 



In analogy with the standard RWA of quantum optics, 
we obtain an approximate Hamiltonian to describe the 
rnth resonance by neglecting all the terms in H m (t) with 

oscillatory time-dependence: H m (t) ps h . This ap- 
proximation is valid as long as the condition 



5 {m \jJ„ 



(19) 



holds (20|. 

Finally, we obtain the time-independent effective Hamil- 
tonian 



;(m) 



! = 1 



i=l 



where the parameters J : 



(m) 



(20) 



and 4 = f [1 - (-l) m J m (^)] denote effective 
anisotropics in the rotating frame. Interestingly, the ef- 
fective Hamiltonian Eq. (|20[) corresponds to an exactly 
solvable model, i.e., it is unitarily equivalent to the XY 
anisotropic spin chain in a transverse field [38M40l | . How- 
ever, in our case, the anisotropics depend upon cither 
the order m of the resonance as well as the driving am- 
plitude gi . Therefore, the driving amplitude of the local 
field now plays the role of a new parameter that influ- 
ences the criticality of the system. 



B. Signatures of criticality in the rotating frame 



Under the RWA, the Hamiltonian and the Flo- 

quet Hamiltonian ti = H(t) — zjj are isospectral oper- 
ators, which implies that the eigenvalues of the effective 
Hamiltonian correspond to the quasienergies. 
As we show in appendix B — similarly to the Ising model 
[l[ , after a Jordan- Wigner transformation, and a discrete 
Fourier transform, the Hamiltonian Eq. ([20| can be writ- 
ten as follows 



;(m) 



]T [(26^ - w*)(4<* + cL k c. k ) - 28^' 

k>Q 

J2(-irA k J m f^-) (clcl k + c-kCk) 

k>0 ^ ' 

E^o?- (21) 



/c>0 



The matrix representation of h m k in the Nambu sub- 
space is given by 



h {m) - 



(22) 

The Hamiltonian Eq. (f2~Tj) can be diagonalized via a Bo- 
goliubov transformation 



fe>0 



E e ^™ 



t. 



(23) 



where 



£k,-. 



2]J (5M - Jcosfc) 2 + (jj m ^ 



(24) 

Furthermore, the quasienergies arc defined (modulus tt) 
by the equation 



,(±) 

"k.m 



= -u> k ± e k . 



~2~' 



(25) 



as defined in Eq. ((9]). The quasienergy gap in the 
fcrmion picture is given by AE k , n 



-(+) J-) 



2e k 
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Therefore, when the gap closes — modulus f2, the effective 
Hamiltonian exhibits a behaviour which resembles the 
dynamics of a critical quantum system. The continuous 
lines in Fig. [1] (b) depict the quasienergy dispersion rela- 
tion for gi ^ 0. We observe that the driving lifts the de- 
generacy giving rise to an anticrossing. Based on the well- 
known results for the time-independent XY model that 
we summarize in appendix B, we find that the system de- 
scribed by the effective Hamiltonian Eq. (|20|) exhibits an 
effective nonequilibrium Ising-like QPT along the criti- 
cal lines |<j( m )| = J, and a nonequilibrium anisotropic 
QPT along the lines where J m (%-) = 0, as long as the 
condition |c)( m )| < J holds. The gapless quasienergy ex- 
citation spectrum for parameters along the critical lines 
is as a direct consequence of coherent destruction of tun- 
neling [23|, H3] , i.e, of the existence of a generalized parity 
symmetry in the extended Hilbcrt space 11® T, where 1Z 
is the Hilbert space of square integrable functions and T 
the space of time-periodic functions [22[ . 

Figure [2] (a) depicts the character of the quasienergy 
excitation spectrum ±Ek, m for parameters in the ferro- 
magnetic phase FMY, along the critical line, and in the 
ferromagnetic phase FMZ, respectively. Figure [21(b) de- 
picts the phase diagram for the nonequilibrium QPT in 
the neighborhood of the 2-photon resonance. The white 
zones in the phase diagram correspond to the effective 
paramagnetic phase and are defined by the inequality 
J < |c5( Tn) | < |^|, for to = 2, where 5^L denotes the 
maximum detuning for which the RWA is still valid. The 
anisotropic transition is characterized by two ferromag- 



netic phases, i.e., for J? 



> J, 



(m) 



the system is in a 



ferromagnetically-ordered phase along z direction FMZ, 
while it is the other way around in the FMY phase. In 
Fig. [5] (c) we plot the the effective asymmetries ji™'' and 
Jy as a function of the driving amplitude g\ in the case 
of a 2-photon resonance. 



III. PHYSICS IN THE LABORATORY FRAME 
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FIG. 2: Nonequilibrium quantum phase transition in the 
driven Ising chain in a time-dependent transverse field g(t) = 
go + gi cos ilt. (a) Depicts the quasienergy dispersion rela- 
tions ±£fc,m for parameters in the ferromagnetic phases FMZ 
and FMY, and along a critical line, (i>) Depicts the phase 
diagram of the nonequilibrium phase transition around the 
m = 2 resonance as a function of the driving amplitude gi 
and the static local field go. The white zones represent the 
paramagnetic phase. Correspondingly, the blue zones repre- 
sent the ferromagnetic phase FMZ and the green zones the 
ferromagnetic phase FMY. (c) Depicts the effective asymme- 



tries in the z direction Ji (blue curve) and y direction J y 
(green curve) as a function of the driving amplitude gi. For 
this plot, we have considered J/Q — 0.01. 



A. Quantum evolution of an initial paramagnetic 
state 



and, correspondingly, the negative-quasienergy Floquct 
mode is given by 



As we previously mentioned, the eigenvalues of the ef- 
fective Hamiltonian correspond to the quasienergies of 
the system. However, the corresponding eigenstates do 
not necessarily correspond to Floquet modes. In order 
to obtain the Floquet modes, one should apply a unitary 
transformation back into the laboratory frame. In so do- 
ing, the positive-quasiencrgy Floquet mode around the 
TO-photon resonance in the Nambu subspace is given by 



cos(4> k: , n ) 

~M<i>k,m) 
e _ <(Wt) _^ t)cofl( j> 



sin(0 feim ) 



(26) 



e _ 4(2am(t) _- t)sin(0fem) 
e . i(2Qm(t)+ ^ t) cos{4>ktm) 



where 



tan(20 fe 



-(-l) m A fc J m (^) 



2<5(^ 



(27) 



(28) 



Now, let us investigate the quantum evolution in the lab- 
oratory frame around the TO-photon resonance when the 
system is initialized in a paramagnetic state of the un- 



6 



driven model with all the spin polarized along the x axis 
IV'm,O) = 0|O_ fcj Ofc>. (29) 

A;>0 

Restricted to the Nambu subspace for a given < k < n, 
such an initial state corresponds to the spinor 4^ m (0) = 
(u* k m (0), vl m (0)) = (0,1), whose quantum evolution is 
given by 



**,m(*) = ?7i m) (t)exp(-i^7t)* feim (0) 



(-) 



k,m > 



(30) 



where = e~* E *>.™*$^(t) denote the Floquet 

states restricted to the Nambu subspace. 



B. The Dynamics of the transverse magnetization 

By using the exact quantum evolution of the initial 
paramagnetic state, we are able to calculate the dy- 
namics of physical observables in the laboratory frame. 
The transverse magnetization M x (t) gives us information 
about the occurrence of a macroscopic polarization of the 
spins along the x axis. Let us consider the expectation 
value M x (t) = jj(tp„ 
resonance, 



t\ J2iLi a i iV'm, t) close to the m-th 



r n dk 

M x (t) = - I \(t)a*(k)* k!m (t) 



7T 



f n dk 

1-2 — sm 2 (e k , m t) sin 2 (2</. fc!in ), (31) 

J0 7T 



where 



a{k) ~2d 9Q ~\0 -1 



(32) 



and n k is the FBdG Hamiltonian Eq. ©. Fig. [3] de- 
picts the dynamics of the magnetization in the thermo- 
dynamic limit calculated using RWA (black curve). In 
particular, Fig. [3] (a) show the dynamics for parameters 
in the nonequilibrium paramagnetic phase and Fig. [3J (b) 
for the ferromagnetic phase FMZ. We observe that in 
the paramagnetic phase the system exhibits a station- 
ary state which corresponds to a polarized state along 
the local field direction. In contrast, for parameters cor- 
responding to the ferromagnetic phase FMZ, the state 
is not totally polarized along this direction. The green 
curve in Fig. [3] depicts the result of exact numerical calcu- 
lation (see appendix C) of the magnetization for a finite 
size system N = 100. 

In this dynamical scenario the connection with criticality 
is not obvious. Rather, signatures of quantum criticality 
in the laboratory frame may appear in the asymptotic 
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FIG. 3: (a) Time evolution of the magnetization M x (t) in the 
thermodynamic limit (black curve) and exact numerical result 
for a finite system consisting of iV = 100 spins (green curve) 
for < t < 200T, where T = 2n/Q is the period of the exter- 
nal driving, (a) Time evolution for parameters correspond- 
ing to the nonequilibrium paramagnetic phase > J), 
for m = 2, where (g 1 /n,g /Q) = (1,0.515). (b) Time evo- 
lution for parameters corresponding to the nonequilibrium 
ferromagnetic phase (FMZ) < J), for m = 2, where 
(gi/Q, ga/Q) = (1,0.505). The insets show the detail of the 
magnetization curves. We have considered J/Q = 0.01. 



dynamics. Let us consider now the time evolution of the 
expectation value of a general observable 



o{t) = (^ m ,t\d\ip m ,t) 



dk 



2tt k ' m 



In general, following the argument established in [331 ] . it 
is possible to show that 0(t) = O pcl '(t) + O tx {t), where 



O peT (t) 



E 



Ae{+,-} 



n dk 

27T 1 



(34) 

is the periodic contribution to the expectation value, 
which correspond to synchronization with the external 
driving. Here we consider A + = — sin[0/ £ m ] and A_ = 
cos[(j)k,m] ■ Correspondingly, 



O tT (t) 



* dk 



3? 



A; ^ e - 2 ^, mt[$ W (0] t Ofc$ (-) W 



(35) 

denotes the transient component, which decays to zero in 
the long-time limit as a consequence of destructive inter- 
ference in time [33J. Therefore, the system tends to syn- 
chronize with the external control in the long-time limit. 
We conclude that in the asymptotic dynamics, the Flo- 
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FIG. 4: Cycle-averaged energy ffm as a function of the 
driving amplitude g±. (a) For parameters in nonequilibrium 
paramagnetic phase (8^ m ' > J), for m = 2, where po/^0.515. 
(b) For parameters inside the ladder of ferromagnetic phases 
(r5 (m) < J), for m = 2, where p /^ = 0.505. The insets depict 
the second derivative of the cycle-average energy as a function 
of the driving amplitude gi. We have considered J/fi = 0.01. 



FIG. 5: Cycle-averaged magnetization Mz as a function of 
the driving amplitude g±. (a) For parameters in the nonequi- 
librium paramagnetic phase (<5' m ' > J), for m = 2, where 
<?o/f2 = 0.515. (b) For parameters inside the ladder of ferro- 
magnetic phases ((5' m ' < J), for m = 2, where <?o/f2 = 505. 
The insets depict the second derivative of the cycle-average 
magnetization as a function of the driving amplitude g\ . The 
parameters are m — 2 and J/Q = 0.01. 



quet modes determine the quantum critical behaviour, as 
we discuss in the next section. 



Cycle-averaged expectation values in Floquet 
eigenstates 



We now define cycle-averaged expectation values of 
physical obervablcs. In the case of a time-dependent 
Hamiltonian H(t), the energy is not conserved. There- 
fore, to describe signatures of the quantum phase tran- 
sition in the laboratory frame we define the aver- 
aged energy Hm in the Floquet state \^m(t)) = 
(±) 



fe >O \n fm(*)> *» 



* dk 



dt 

T Jo 2tt 
T dt r dk 
TJo 2^ 



4 ± 2 l + (Cri(*)) t C^W 



d 



,(±) 



(9/ 



' k,m \ 



(36) 



By using the analytic expressions for the Floquet modes 
Eqs. (j2"o) and (f2"T)) we obtain the expression 



2tt 



— cos(20 fcjm ) 



(37) 



On exact resonance 5( m ) = 0, we obtain an analytical 
expression for the cycle averaged energy 



7T 



1 



4gi 



(38) 



where is the complete elliptic integral of the second 
kind (see appendix B). This result confirms our predic- 
tion based on the description of the system in the ro- 
tating frame (see Fig. [2]). The cycle-averaged energy ex- 
hibits singularities at the zeros of the Bessel function, 
i.e., J m (-$f) = 0. This is a clear signature of a criti- 
cal nonequilibrium behaviour. Finally, we calculate the 
cycle-averaged transverse magnetization in the Floquet 
mode (t))- By considering the extension of the 

Hellmann-Feynman theorem for Floquet theory [U [22J , 

we can compute the cycle-averaged magnetization 
in terms of derivatives of the quasienergies as follows 



Mi±> = - 



T dt r dk , 



T . 
T dk d, 



2tt 

'k,m 



dgo 



2?r dgo 



dk 

— cos(20 fc;m ). 

7T 



(39) 



Figures 0] and [5] show the cycle-averaged expectation 
values of observables. Similarly to the undriven case, the 
system exhibit nonanalyticities in the second derivative 
of observables — as can be seen in the corresponding in- 
sets, thus resembling a continuous phase transition. Such 



8 



nonanalyticities arise as a consequence of CDT [23J, \2A 
and therefore, from the gapless quasienergy spectrum. 



IV. CONCLUSIONS 



We have investigated the nonequilibrium critical be- 
haviour in a driven 1-dimensional transverse Ising model. 
We find that the fundamental signature of critical be- 
haviour is the existence of a gapless quasienergy spec- 
trum, which is a direct consequence of CDT. The role of 
coherent destruction of tunneling in nonequilibrium QPT 
has been explored either theoretically as experimentally 
in the context of driven superfluidity [HI, [l4j]. In this 
paper, we show that CDT induces a critical behaviour 
which resembles a second order QPT. In particular, the 
symmetry which is broken corresponds to a generalized 
parity in an extended Hilbcrt space where TZ is the 

Hilbcrt space of square integrable functions and T the 
space of time-periodic functions. In this nonadiabatic 
scenario, the short time dynamics is governed by tran- 
sient effects that tend to zero in the asymptotic limit as 
a consequence of destructive interference [33[. The long 
time dynamic is governed by the Floquet modes, there- 
fore, the nonequilibrium quantum phases correspond to 
states which arc synchronized with the external control. 

Possible experimental implementations of our model 
could be achieved using cold atoms [U, H, @, and in where 
fully C-labelled sodium butyrate using liquid state nu- 
clear magnetic resonance (NMR) [42|, EH ■ 
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the Jordan- Wigner representation of the Pauli matrices 



< 1 I'',':,- 



J-1 



3-1 



^■={c]+c ] )\{(l-24 



CI J 



(Al) 



Under this representation of the angular momentum al- 
gebra, the parity operator Eq. ([3]) acquires the form 



N 



n = II(i-2ctcj 

3=1 



(A2) 



Using this parity operator we are able to define projectors 
on the subspaces with even (+) and odd (— ) number of 
fermionic quasiparticles as follows 



p± = -(i±n). 



(A3) 



The projectors V± satisfy the usual properties of or- 
thogonal projection operators such as V+ + V- = 1, 
V+P- = V-V- = 0, and (V±) 2 = V±. Using this prop- 
erties and the algebra of fermionic operators is possible 
to show that the Hamiltonian Eq. ([1]) can be decomposed 
as follows 



H(t) = V + H (+) (t)V++V-H i - ) (t)V- 



(A4) 



N 



N 



H&(t) = - 5 (t)^(l-24c i )-J^(cJ-c i )(ct fl +c i+1 ). 

i=l i=l 

(A5) 

To perform the splitting we have defined antiperiodic 
boundary conditions in the even (+) subspace cn+i = 
—C\ and periodic boundary conditions in the odd (— ) 
subspace cyv+i = c\. In this paper we focus on the pro- 
jection in the even subspace. Translational invari- 
ance suggest to use the discrete Fourier transform 



-in/4 

\ * _ ikn 



c k e 



(A6) 



which is compatible with the antiperiodic boundary con- 
ditions when k £ {±f ,±^,...,± (jV ~ 1)7r }. The discrete 
Fourier transform maps into Hamiltonian Eq. (jj). 



Appendix A: Description of the periodically-driven 
Ising model for finite size 



Appendix B: The quantum phase transition in the 
anisotropic XY spin chain in a transverse field 



In this section we introduce the fundamental tools used 
in the solution of the Ising model following the methods 
and the notation of reference J25| . 

Let us consider the Hamiltonian Eq. (fTJ) in the case of 
even number of lattice sites N. For convenience, we use 



Here we consider the critical behaviour in a model de- 
scribed by the Hamiltonian 



iY 



H = -^<-^(J^^ +1 + J y a,faf +1 ),(Bl) 
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which is unitarily equivalent to the Hamiltonian of an 
anisotropic XY spin chain in a transverse field (38l - l40j . 

Similarly to the Ising model, after Jordan- Wigner 
transformation, and a discrete Fourier transform, the 
Hamiltonian Eq. (|B1[) in the even subspace (the sub- 
space with a even number of fermionic quasiparticles) 
acquires the form 



where E[z] is the complete elliptic integral of the second 
kind [13] . The scaled ground state energy exhibits a non- 
analyticity of the second derivative at the critical line of 
the anisotropic transition 7 = which is a generic 
characteristic of a second-order QPT 

Appendix C: Numerical calculation of the 
expectation values 



H = Y^{ 2 ( h -( J * + J y) cosk )( c l c k + c -k c -k)- 2 h 

k>0 

^2(J Z - J„)smfc(4cl fc +c_ fc c fc ). (B2) 



fc>0 



The diagonalization of this Hamiltonian is completed af- 
ter a Bogoliubov transformation, 



(B3) 



where 



E k = 2y (h - ( J z + J y ) cos k) + {{J z — J y ) sinfc) . 

(B4) 

The system exhibits an Ising-like QPT along the lines 
\h\ = J z + J y and an anisotropic QPT along the line 
Jz = Jy> providing that \h\ < J z + J y . The anisotropic 
transition is characterized by two ferromagnetic phases, 
i.e., for J z > J y the system is in a ferromagnctically- 
ordered phase along z direction FMZ, while it is the other 
way around in the FMY phase. Similarly to reference 
[381 ] , we consider a reparametrization of the asymmetries 

Jz = ^(1 +7)) 



■ In 



-(1-7), 



(B5) 



By using the BCS ansatz Eq. ([S]) we can solve the 
Schrodinger equation for the Ising model in terms of 
the solution of the Schrodinger equation for an effec- 
tive two-level system described by the BdG Hamilto- 
nian Eq. (IH]), which is parametrized by the quasimo- 

}. In the numer- 



mcntum k € {±77, ±4£, • • • , ± 



N 



ical calculation we assume that the system is prepared 
initially in the unoccupied state |0_fe,0fc), which implies 
that #£(0) = (it*(0),w*(0)) = (0,1). After numerical 
integration of the dynamical BdG Eq. (J7J, we find the 
spinor ^fc(t). To calculate the scaled expectation value 

of the transverse magnetization M x (t) = (^YliLi <T \ 
for a given system size N we use the formula 



M x (t) 



N 



X)(**(t))V(*)]**(*)- 



(CI) 



fc>0 



In the last expression we have used the definition of a z (k) 
given in Eq. (|32|) . For example, to calculate the dynamics 
of the system for N = 4, we perform the numerical inte- 
gration of Eq. (0 for k <E {f , ^}. Based on the solution 
of this equation we find the solution of the Schrodinger 
equation for Hamiltonian Eq. ([T]) using the BCS ansatz 



where 7 is dimensionless parameter characterizing the 
degree of anisotropy in the zy-plane. Under this 
reparametrization, the Ising-like critical lines correpond 
to \h\ = J, and the anisotropic transition occurs at 7 = 0, 
as long as \h\ < J. 

Interestingly, in the absence of a transverse field, i.e., 
for h = 0, the scaled ground-state energy can be written 
in the thermodynamic limit as 



Er 



lim 



E k 



= -j 



N 
dk 



7 2 ) sin 2 k 



T^ 1 

7T 



(B6) 



|V, t) =|^ f ,t)®|VaL,t), 



where 



and 



(C2) 



\^,t) = (u f (t)|l_ } ,l f ) +v i (t)\0- l ,0 f }) , (C3) 



\ij>^t) = (t»to(t)|l_a P L,la Jt )+i;a J t(t)IO_a t ,Oa |t ) 



(C4) 
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